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FINITE -ELEMENT COMPUTER PROGRAM FOR AXISYMMETRIC LOADING 
SITUATIONS WHERE COMPONENTS MAY HAVE A 
RELATIVE INTERFERENCE FIT 


by Christopher M. Taylor* 

Lewis Research Center 

SUMMARY 

This finite -element computer program permits the analysis of axisymmetrically 
loaded components having a relative interference fit. External, inertial, and thermal 
loads may be applied in addition to forces arising from an interference condition. The 
components under investigation may comprise different homogeneous isotropic materials 
whose properties may be temperature dependent. The theoretical background to the anal- 
ysis of the interference condition is described, and the computer program structure out- 
lined briefly. Particular attention is paid to the program input requirements and the out- 
put. The output comprises the input information, nodal distortions, element stresses, 
and details of the interference contact conditions. The running of multiple cases, with 
one or a number of different geometries, is possible. An example to illustrate the data 
input and program output is presented. 


INTRODUCTION 

The ability of the engineer to undertake the analysis of deformations and stresses 
occurring in the structures and components with which he deals has changed remarkably 
in the last 15 years. The development of finite -element approximation techniques, cou- 
pled with the widespread availability of digital computers, has paved the way for detailed 
examination of the behavior of loaded structures and machine parts of quite complex 
geometry. Numerous computer programs have been developed by workers in industry, 
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research organizations, and institutions of higher education using the finite -element 
technique. Some of the more recently developed programs enable the study of situations 
in which three-dimensional loading pertains. Work to improve the accuracy of predic- 
tion in such situations, whilst most efficiently employing the facilities of the digital com- 
puter, progresses currently. 

The program described here restricts itself to simple elements and axisymmetric 
loading situations. Unlike most finite -element programs, however, it can undertake a 
direct analysis of components that have relative interference. It is a frequent engineer- 
ing practice to mate or prestress parts using shrink or press -fitting techniques. The 
analysis of such situations has, in general, been limited to axially invariant axisymmet- 
ric situations using, for example, the solution for a cylinder subjected to uniform inter- 
nal or external pressure due to Lame (ref. 1). 

Wilson and Parsons (ref. 1) have described the approach adopted here to enable the so- 
lution of problems in which statically indeterminate conditions occur as a result of inter- 
ference fits. They term the technique "differential displacements. " It enables one node 
in a finite -element mesh to be displaced by a specified amount relative to another without 
having to know how the two nodes move in relation to some fixed coordinate system. 

Since the finite elementry used in the computer program developed is quite conven- 
tional, the details are not dwelt upon. A number of texts (e.g. , ref. 2) may be consulted 
for background information. The differential displacements approach and its incorpora- 
tion into solution schemes for systems of equations is described. The computer program 
is presented as an appendix. As with many finite- element programs, it is long and its 
details are not discussed; its broad structure is, however, described. Fuller attention 
is paid to the input data for the program and the output to enable its easy use. An ex- 
ample of the inner ring of a roller bearing shrunk or press fitted onto a shaft is pre- 
sented. The purpose of this example is not to highlight the technical aspects of the sit- 
uation, but rather to give examples of the input and output. 

The computer program presented uses some coding developed by Dr. B. Parsons, 
Reader in Mechanical Engineering, University of Leeds, England, whose encouragement 
and support are gratefully acknowledged. 


SOME ASPECTS OF THE ANALYSIS 

The finite -element types, designated 1, 2, and 0, that have been incorporated into 
the program are shown in figure 1. The geometrical restrictions of each element form 
are indicated. A linear displacement model within the elements was assumed, enabling 
the stiffness characteristics in matrix form to be determined by analytic integration of 
the strain energy integrals. The unknowns are normally node displacements in the 
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radial and axial directions. Where a node displacement or a coordinate direction is pre- 
scribed, an unknown force results. The general equation enabling nodal distortions to be 

determined is 


[K]{ U } = { F } W 

where [K]is the overall stiffness matrix of the structure, {U} is the vector of nodal 
displacements, and { F} is the vector of nodal forces. The solution of the system of 
equations (1) may be achieved by many schemes. Direct methods of solution are now al- 
most universally employed for solid mechanics problems, the method used here being 
that of Crout (ref. 3). 

If it is known that two nodes in a finite-element structure have an effective pre- 
scribed movement relative to each other (obviously an interference), it is cumbersome 
with conventional finite elementry to analyze the situation. An influence coefficient ap- 
proach can be adopted (e.g., ref. 4), but to use the differential displacement approach is 
more direct and less specific. 

For example, consider, for convenience, consecutively numbered nodes 1, 2, 3, 
and 4 where it is known that nodes 2 and 3 have an interference { 6} . With due regard 
for sign we may write, 


{ Ug } = {U 2 } + {6} 

{F 2 } = {F’} - {F 6 } 
{F 3 } = {F’} + {F 6 } 


( 2 ) 


(3) 


where {F 2 } and { Fg } are node forces due to external, inertial, and thermal loadings 
and {F g } is the interference force due to the differential displacement {6} and is 
unknown. 

Thus the equation of the form of (1) for the nodes may be adapted to 
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Equation (4) can be reduced to give 
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Having made this reduction, the set of equations (5) can now be solved for the unknown 
displacements, all other quantities being calculable. After the determination of the dis- 
placements, the set of equations (4) may be used to calculate the nodal forces due to in- 
terference { F^} . 

The differential displacements technique, already described in terms of four con- 
secutively numbered nodes, two of which have relative interference, is amenable to pro- 
gramming on a digital computer for large systems of equations. Nodes having inter- 
ference need not be numbered consecutively. The example presented in the report shows 
how the geometry of interferring nodes is described. 


THE COMPUTER PROGRAM 

The program listing is presented in appendix A. The program is written in Fortran 
and has been developed on a Univac 1100 machine with a Fortran V compiler. In appen- 
dix B the structure of the program is briefly described. To facilitate the use of the pro- 
gram, detailed attention will now be given to the required input and output. 


Input Data 

The required card input will be listed in sequence with comments on limitations. 
Some general remarks will first be made. Where detailed elucidation on some aspect of 
the input data is needed, this appears in appendix C. 

The coordinate system adopted is shown in figure 2. For elements of types 1 and 0 
the axial node positions (which are relative to node i) should be in accord with the sign 
convention of figure 2 (i.e. , values to the left of node i being negative). 

Any consistent set of units may be employed. The program output does not detail 
any units, but those actually employed could be indicated in the title. 

The program arrangement permits multiple cases on one run with the same geom- 
etry but varying operating conditions, for example, loading or speed. In addition 
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different geometries may be examined on the same run, though this is a less likely 
requirement. 

The input requirements of finite -element computer programs are usually demanding 
in terms of punching of data. To simplify this aspect with the present program, it has 
been assumed that the most likely situation is that a single material (termed the "base” 
material) will be involved or will dominate. Other materials (termed "special"), un- 
limited in number, may be used in addition to the base material. However, when many 
special materials are used, the preparation of data cards will be disproportionately 
more time consuming than for a single material case. 

The program includes a consideration of external, inertial, and thermal loadings. 
Material properties may be a function of temperature. The temperature distribution is 
an input requirement. 

The limitations on the magnitude of certain variables indicated in the required input 
relate to the dimensioning of arrays. Alterations may be made as indicated in appen- 
dix C, where the storage requirement needed for operation on a Univac 1100 computer 
is also given. 

It should also be noted that axisymmetric loading problems in which no interference 
fit occurs can also be treated with the program. The program variable for a parameter, 
on read -in, is indicated in brackets after its verbal description. 

CARD A format (12) 

12 number of geometries considered (IC) 

The remaining cards are repeated for each geometry. 

CARD B A title of up to 80 characters, commencing with a blank 

CARD C format (12) 


12 number of different cases for geometry under consideration (NOC) 

CARD D format (614) 

14 number of elements (M) 

14 number of nodes Cf200) (N) 

14 maximum nodal number difference in any element Cf>19) (NW) 

(See appendix C . ) 

14 number of special material elements (NE) 

14 number of nodes constrained 0^50) (NCON) 

14 number of nodes with external forces applied (NFOR) 
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CARD E format (5E12.4) 

E12.4 Poisson's ratio for base material (PROP (1)) 

E12.4 Young’s modulus for base material (PROP (2)) 

E12.4 coefficient of expansion for base material (PROP (3)) 

E12.4 density of base material (PROP (4)) 

E12.4 reference temperature (<£o) (PROP (5)) 


If any of the first four properties read in on this card are to be functions of the fifth, 
temperature, then any negative number should be placed in the appropriate format posi- 
tion. The reference temperature is that temperature the elements will be assured to 
take if the element temperature of CARD H is left blank. This reference temperature 
may be zero, which implies no relative thermal expansion for the element; that is, if the 
geometry read in is to suffer no thermal distortions, the reference temperature should 
be zero and apply to every element. 

CARD F format (F10. 0) 

F10. 0 rotational speed (rpm) (W) 

CARDS G,H 

CARD G format (614) 


14 

element number 

(L) 

14 

type of element 

(NTYPE) 

14 

number of node i 

(NTRI (L, 1)) 

14 

number of node j 

(NTRI (L, 2)) 

14 

number of node k 

(NTRI (L, 3)) 

14 

number of node L (0 for a triangle) 

(NTRI (L, 4)) 
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CARD H format (6F1 0.0) 


Element type 



1 

2 

0 


F10.0 

r i 

r i 

r i 

(TRI (L,l)) 

F10.0 

r j 

” r j 

r j 

(TRI (L, 2)) 

F10. 0 

b i 

r k 

b i 

(TRI (L, 3)) 

F10.0 


r L 

b k 

(TRI (L, 4)) 

F10.0 

b L 

b 

0 

(TRI (L, 5)) 

F10.0 

Element temperature (<£o) 

(TRI (L, 6)) 


Figure 1 may be referred to for the meaning of symbols on cards G and H. If the ele- 
ment temperature is to equal the reference temperature of card E, a blank space may be 
left, otherwise the temperature of all elements must be indicated here. Note that if any 
temperatures are to be zero then the reference temperature must also be zero. 

Cards G and H are repeated, consecutively, for each element. For type 1 rectan- 
gular elements only, element numbering need not be consecutive. In this case inter- 
mediate elements will be generated automatically. The automatic generation defines 
node numbers of elements by increasing the previous element values by unity; linear 
interpolation is used to determine r., r^ and element temperature; and b., and b L 
remain fixed. The automatic generation may be applied in either coordinate direction. 

CARD I format (14, 4E12.4) Only read if number of special material elements (NE) 


of card D is greater than zero 

14 number of special material element (L) 

E12.4 Poisson's ratio for special material (TRI(L,7)) 

E12.4 Young's modulus for special material (TRI(L, 8)) 

E12.4 coefficient of expansion for special material (TRI(L,9)) 

E12.4 density of special element (TRI(L,10)) 


The properties should have a value appropriate to the element temperature of card E or 
card H. Card I should be repeated "NE times" corresponding to the number of special 
material elements. 
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CARDS J,K (For cases when properties of base material are temperature dependent) 
CARD J format (12) 

12 number of temperatures at which property value given (>10) (IPROP (I)) 
CARD K format (F10. 0,E12. 4) 

F10.0 temperature (U (J)) 

E12.4 value of property at this temperature (U (K)) 

Card K is repeated for the number of temperature values indicated on card J, with tem- 
perature values increasing. 

Cards J and K are repeated for any property of card E that was read in as a negative 
number. The cards should be presented for properties in the same sequence as indicated 
for card E (i. e. , Poisson’s ratio, Young's modulus, coefficient of expansion, density) 
and omitted entirely if a property is not temperature dependent. The range of temper- 
atures on card K must encompass the full range of element temperatures. 


CARD L format (314, 2F 10. 4) 

14 number of node (NUV (1, 1)) 

14 = 1 if node movement constrained radially, otherwise 0 (NUV (I, 2)) 

14 = 1 if node movement constrained axially, otherwise 0 (NUV (I, 3)) 

F10.4 amount of radial movement if constrained (ANUV (1, 1)) 

F10.4 amount of axial movement if constrained (ANUV (1,2)) 


Normally the constrained movement will have zero value. Nodes on the axis of revolu- 
tion must be restrained to zero radial movement. Cases of axial symmetry can be ac- 
commodated considering only half the problem by using appropriate axial constraint. 
Card L should be read in for NCON nodes C}>50), as indicated on card D, and omitted if 
no nodes are constrained. 

CARD M format (I4,2E12.4) 


14 

number of node 

(L) 

E12.4 

applied radial force 

(PROP (1)) 

E12.4 

applied axial force 

(PROP (2)) 


Card M should be repeated for NFOR nodes as indicated on card D and omitted if no 
nodes have applied forces. An applied distributed loading must be appropriately divided 
between nodes, the force at each being that for 2j t radians (see appendix C. ) 
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CARD N format (12) 

12 = 1 if there is any interference, otherwise 0 (IDD) 

If card N reads in the value 0, the following course is adopted: if the number of 
cases for the geometry under consideration is 1 (NOC of card C), then the data ter- 
minate, unless another geometry is to be considered, in which case the card sequence 
returns to card B and begins again. If there are (NOC - 1) cases for the same geometry 
still to be dealt with, then the following cards should be read for each of these cases: 
CARD F, CARD M, and CARD N. This sequence allows applied speed and load to be 
varied. If card N should now take the value 1, then the case should proceed as below. 

If card N reads in the value 1 (an interference fit), the following card sequence 


should be adopted: 

CARD O format (15) 

15 number of node pairs with relative interference (NODD) 

CARD P format (IX, 16, 5X, 16, 5X, 16, Ell. 4, 5X, 16, Ell. 4) 

16 lower numbered node of a pair in interference (J) 

16 higher numbered node of a pair in interference (K) 

16 = 1 if radial interference, otherwise 0 (HR) 

Ell. 4 amount of any radial interference (DR) 

16 = 1 if any axial interference, otherwise 0 (12) 

Ell. 4 amount of any axial interference (D2) 


Card P should be repeated for the NODD pairs of nodes indicated on card 0. The 
sign of interference, radial or axial, should be in accord with the sign convention of fig- 
ure 2, working from the lower to the higher numbered node. (See the example of ap- 
pendix D. ) It will be convenient in interpreting the output to have all the higher num- 
bered nodes of pairs in radial interference "nearer M the axis of revolution. This 
scheme has been adopted in the development of the program and is recommended. 

CARD Q format (14) 

14 = 1 if radial interference pressures are to be calculated (IND) 

If card Q reads in 1, then the pressure distribution at a fixed radial location where 
interference occurs will be determined. Card Q may only be read in with a value of 1 if 
all the nodes having interference are at the same radius and are axially adjacent. If 
this card reads in 0, go to card T. 
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CARD R format (14, E12.4) 

14 number of node pairs in interference Cj-50) (NP) 

E12.4 constant radius of interference fit (RAD) 

CARD S format (I4,E12.4) 

14 number of the node in pair nearer the axis of revolution (IP (I)) 

E12. 4 axial distance to next node in sequence AL (I) 


Read in the node numbers nearer the axis of revolution in sequence, moving axially in 
either direction. The axial distance between nodes should be kept positive. The pres- 
sure distribution is the same for both components involved in the interference. The 
value of axial distance to the next node for the last node is immaterial. 

Card S should be repeated for the number of node pairs (NP) indicated in card R, 

CARD T format (12) 

12 = 0 if no axial slip to be allowed automatically at nodes with radial 

interference, otherwise 1 (1SLP) 

This card may only read in a value of 1 if it is permissible to do so with card Q. If card 
T reads in a value of 0, the following card sequence should be adopted for remaining 
cases with the same geometry: CARD F, CARDS L,M,N, et seq. Alternatively, if 
there are no more cases for the particular geometry, the data sequence terminates or 
returns to card B for a new geometry. 

If card T reads in a value of 1, the following card should be inputted: 

CARD U format (F5. 0) 

F5. 0 coefficient of friction (COF) 

See appendix C for an explanation of the slip calculations. After this card, the data 
sequence should be the same as if card T had read in 0 (see above). 

Although it is very convenient and saving in time to run multiple cases, particularly 
with the same geometry, it will be seen that the data is accordingly complicated. It is 
therefore recommended that new users of the program adopt single case runs initially. 
The data card series will then be as in table I. 

Output 

The program output is virtually self-explanatory. An example is to be found in ap- 
pendix D. As well as outputting calculated information, all the input data are also pre- 
sented. Thus, the first part of the output gives the data fed to the program on cards B, 
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D through M, and P, as appropriate. Calculations to set element temperatures and 
physical properties are carried out on the basis of the above input but details are not out- 
putted. In the case of multiple runs with the same geometry only the parameters which 
can be changed are re-presented. 

Results calculated by the program and presented in the output are 

(a) Nodal deformations 

(b) Element stresses at the midpoint of the element. Stresses given are radial, tan- 
gential, axial, and shear together with the coordinates of the element midpoint with re- 
spect to node i. Positive stresses are tensile. 

(c) If there is an interference fit, the forces due to this fit radially and axially acting 
on the higher numbered of each node in the pairs involved. The sign of these forces is 
opposite to the sign convention of figure 2. Thus in the example of appendix D, where all 
the higher numbered nodes are nearer the axis of revolution, the positive radial forces 
indicate that the two components are still entirely in interference. These forces are per 
radian unlike the input forces, which are for an extent of 2jt radians. 

(d) If the option is so chosen, the pressure distribution at the position of radial in- 
terference. Note that although the interference force output may indicate that all nodes 
of an interference remain in contact, small negative interference pressures may some- 
times occur. This is because of the assumption that the pressure varies linearly be- 
tween nodes. 

(e) If slip between nodes in radial interference is to be accommodated automatically, 
the coefficient of friction and details of which nodes are allowed to slip on each cycle . 

On completion of slip the final nodal deformations and element stresses are given. 


CONCLUDING REMARKS 

A finite -element computer program that can directly analyze the behavior of axisym- 
metrically loaded components having an interference fit has been described. External, 
inertial, and thermal loads can be accommodated as well as loading due to the interfer- 
ence fit. In addition components constructed of a number of materials, the physical 
properties of which may be temperature dependent, can be dealt with. The program code 
is presented together with an example data input and program output. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 15, 1977, 

505-04 
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APPENDIX A 


THE COMPUTER PROGRAM LISTING 


C AXISYMMETRIC FINITE ELEMENT PROGRAM WITH VARIABLE 

TEMPERATURE EFFECTS AND INERTIAL LOADING (CMT) 
DIFFERENTIAL DISPLACEMENTS INCLUDED 

DIMENSION ALLK(400, 41), TRI( 20 0,10) ,NT”I (200 ,4 ) ,NCDIS ( 4 ! 
IP ROP ( 5 ) , U ( 4 0 0 ) , NUV(5C,3), A NU V ( 50 , 2 ) , TU ( 400 ) 

2, DEL (400 ) ,DD (400>,NDD (400) ,MOVE (400) , ALLK1 (400 ,41 ) 
3,ISTR(50,4),STR(50,2),N0DIS1(400),IP(E0),AL(5N) 
4,IPR0P(5) 

INTEGER REPT 

COMMON/BLOCK1/TEMP,PR,YM,ALPHA,PHO 
C0MN0N/4L0CK2/M0RD,NW221,NW22,NW 
IC = 0 

READ(E,1G0)ICASF 
1 T C- IC + 1 

WRI TE (6 , 1 c 0 ) 

READ ( 5,101) 

WRITE (6,101) 

WRITE (6 , 1 52 ) 

READ(5,1C0)NCC 
IC1=1 

READ ( 5 , 102 )M, N,Nw ,NE ,NCON ,NFO p 
WRITE (6,154 )M ,N , N W , N E ,NCON,NFOR 
WRITE ( 6 , 1 r 2 ) 

M OP D = ?*N 
NW22 = 2 * N W ♦ 2 
NW221 = NW22+1 
READ( C ,1C3) ( P R0 D ( I ) ,1:1,5) 

READ( C ,1U4)W 

WRITE (6,1 55 ) ( PROP (I ) , 1=1 , 5) , « 

WRITE <6,l c 2) 

W = W*0. 10471976 
CALL SET2 1 ( TR 1 , 2u0 , 10 ) 

CALL SET23(NTRI ,200,4 ) 

WRI TE ( 6 , 1 56 ) 

WRITE (6 , 152 ) 

L ST :C 

DO 2 I = 1 , M 

READ(5,102)L,NTYPE,(NTRI(L,J),J=1,4) 

READ ( 5,105) (TRI (L,J) , J=1 ,6) 

L 1=LST>1 

I F ( L 1 .EQ.DGCT03Q 
L 2 : L - 1 
L3=L-LST 

ZZ= (TRI (L ,1 )-TRI(LST, 1) )/L3 
YY:(TRI(L,6)-TRI(LST,6))/L3 
D 03 1 J1=L1,L2 
N TP I ( J1 ,1 )=NTRI (lST , 1 )*1 
NTRI(J1,2)=NTRULST,2) + 1 
NTRI(J1,3)=NTRI(LST,3>+1 
NTRI(J1,4)=NTRI(LST,4>+1 
TRI ( J l T 3 ) -TRI (LST, 3) 


MATERIALS , 
-0) , FW ( 2C0 ) , 
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IFCABS(Q) • L T • 1 .E-05 ) GOT 08 
IF(C.GT.U.O)GOT^ 9 
Ml r N 1 + 39 
M2 = Ml + 1 

M3 ~ N 1 - 1 

TR I ( I t K ) -U ( M 1 ) +(U(K2) -U(M1))*(TRI(I,6) -U(M'))/(U(M) 
G0T03 

9 CONTINUE 

8 CONTINUE 

6 CONTINUE 

36 CONTINUE 

CALL ^ E T Z 1 (ALLK , 4 C 0 , 4 1 ) 

CALL SE T2 ? ( F W ,200) 

CALL SETZ2(U , 4L0 ) 

CALL SETZ* (N0DIS,400) 

DO 11 L = 1 , M 
TEMP = TRI (L ,6) 

PR = TR I ( L , 7 ) 

YM = TRI ( L f 8 ) 

ALPHA = TRI ( L , 9 ) 
p H0 = TRKLvlG) 

CALL PECT(L,FW,NTRI,ALLK,U,TRT) 

11 CONTINUE 
PEPT- 1 

DO 12 1 = 1 , N 
J = 2*1-1 
J1 = J +1 

ALLK ( J,NW?21 ) = U(J) + FW(I)*W**2 
4 LL K ( J1 ,NW221 ) = U ( J1 ) 

TU( J)=U( J) 

TUC J1 >=U( J1 ) 

U ( J ) = 0.0 

12 U(Jl) : 0 .0 
WRITE (6,176) 

WRITE (6,165) 

IF (NCON*EC.0)GOT08G 
D 0 1 3 I = 1,NC0N 

P EAD(5,1C8) (NUV(I,J),J=1, 7 ),(ANUV(1,J),J=1,2) 

WRITE (6,166) (NUV(I,J) , J = 1 , 7 ) , (ANUV(I , J ) , J = 1 f 2) 

NOD = 2*NU V (1,1) -1 
U(NOD) = ANU V (1,1) 

IFC NUVCI f ?) .GT.NODIS(NOD) >NOOIS (NOD) = NUV (1 , 2 ) 

NOD = NOD + 1 
U (NOD ) = ANU V (I , 2 ) 

13 NODIS(NOD) = NUVCI, 3) 

80 CONTINUE 
WRITE (6,152) 

14 WRITE (6,177) 

WRITE (6,167) 

IF(NF0P.E0.0)G0T081 
D 01 5 1 = 1 , NF OR 

READ(5,1Q6)L, (PPOP(J) , J = 1 , 2 ) 

WRITE ( 6 , 1 68 ) L , (PROP ( J ) , J=1 , 2 ) 

K = 2* L - 1 
KK = K +1 

ALLK C K ,N*221 ) =ALLK CK ,NW221 ) +PROP (1>*0 • 159154^43 

15 ALLK (KK,NW221) = A LL K ( K K , N W 2? 1 ) + PR0P<2>* D. 159154943 

81 CONTINUE 
I DD =G 


“U (M 
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PEAD ( 5 , 1 GO ) IGD 
D 02 0 I=1,M0R0 
0 C2C J- 1 , NW22 1 
20 ALLK 1 (I v j)=ALLK(l t J) 

D 03 7 1=1,400 
37 NCDIS1 (I ) =N00IS(1) 

K TrMO^D 
k CT-C 
T C Y - G 
I WP=0 
0 C T 0 3 55 

39 D C4 C I = 1 # 'J0G 

NODIS (I ) -NODI Si < I ) 

4 C U ( I ) = n 

DC41 I=1,M0RD 
DC41 J=1,NW221 

41 ALLK(I,J)=ALLK1(I,J> 

CALL DIFDIS(KT,ALLK,NODIS,MOVF,N0D,DD,DEL,ISTR,ST r ’,NOnD,KC7) 
GCT047 

33 IF(IDD.GT.O) 

1C ALL DTFDIS (KT f ALLK ,NODIS KDD »D0 f DEL,ISTR , ST^tNOOOtKCT) 

42 CALL BCROUT (ALLK,U, MCRD ,REPT, Kk ,N0uIS ) 

IF( IDP.GT .0 ) GO TO 22 

WRITE ( 6 9 1 50 ) 

WRITE(6,178) 

WRITE (6,169) 

0016 1 = 1 ,* 

* 1 = 2 *T -1 
** 2 = K 1 + 1 

16 WRITE (6,168 > I ,L( hi > ,U(K2) 

WRITE ( 6 , 1 c 2 ) 

WRITE (6,179) 

WRITE (6,170) 

TOTA-P.O 
ST-P.n 
D 01 7 L = 1 , M 
TEMP- TRI ( L, 6 ) 

PR-TRT (L ,7) 

YM-TRI(L,S) 

A LPHA - TR I (L ,9 ) 

17 CALL PECTST (L ,ST,TOTA ,U ,TRI ,NTRI) 

ST-ST / TOT 

WRITE (6,171 ) ST 

IF ( IC1 •EQ*NOC )GOTO 18 

IC1 = IC1 + 1 

READ( 5,104) W 

WRITE < 6 , 1 72 ) W 

W=W*0 . 10471976 

D 01 9 1 = 1 , N 

J=2*I-1 

J1=J+1 

4LLK(J,NW221)=FW(I)*W**?+TU(J> 

19 ALLK(J1,NW221 )=TU(J1) 

P tPT - -1 
WRITE(6,152) 

GOTO 14 
22 KT=MO?D 

D023II=1 , K T 
I-KT-II +1 
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IF CNDDCI >.GE.C) G0T024 
I1=-NDD(I) 

II- MOVE (II) 

0 EL ( I ) -D£ L ( I ) +0(11) 

GCT02 7 

24 I 1 = MO VE ( I ) 

0 EL ( I ) -U ( 1 1 ) 

23 CONTINUE 

46 CONTINUE 

IFCICY.NE .0)G0T045 
WRITE (6,150) 

WRITE(6,178> 

WRITE (6,169) 

D 02 5 1 = 1 , N 
NOD -2*1-1 

25 WRITE (6,168) I, DEL (NOD) ,DEL(N^D + 1 ) 

WRITE(6,152) 

WRITE (6,179) 

WRITE (6,170 ) 

TCTA-O.D 
S T= G • 0 
0029 L = 1 , M 
TE M P- TRI ( L ,6 ) 

PR-TRI(L,7) 

YN=TRI(L,R) 

ALPHA=TRI (L ,9 ) 

29 CALL RECTST(L,ST,T0TA,DEL,TRI,NTRI) 

ST= ST/TOT A 
WRITE (6,171 ) ST 
45 CONTINUE 

IF ( I W R • N E •0)G0TO43 
CALL SET22(DD ,40U) 

0C26 1 = 1, KT 
MC1=I +1-NW22 
W C7=I+NW22-1 
IF (MC1.GT.1 ) J1-HC1 
IF (MC1.LE.1) Jl=l 
IF (MC2.LT.M0RD > J2=MC2 
IF ( M C ? ♦ G E •M0RL)J2-M0RD 
DC27 Jr Ji , J2 

IF(I*LT*J) r )D( I)=DD(I)+ALLK1 (I,J-I+I)*DEL(J) 

2 7 IF(I.GE.J)DD(I)=l‘DCI)+ALLK1 U , I -J + 1 ) *DEL ( J) 

26 CONTINUF 
WRITE (6,152 ) 

WRITE (6,173 ) 

WRITE (6,152) 

D 02 6 I - 1 , N 
MCI =2*1-1 

M C2 =M C 1 + 1 

IF (NOD (MCI ) .GE.0.AND.NDD(MC2) .GE .0)G0T028 
IF (NCD (MCI ) • L T •Q)DD( MCI ) = A L L K 1 (MCI, NW 221 ) - D D ( MCI ) 
IF (NDD(MC1 ) .GE .U)DD ( MC 1 ) =0 . 

IF (NOD(M f2 ). LT. 0) DD (MC2)=ALLK1(MC2,NW 22 1)-DD(MC 2) 
IF ( NDD (MC2) .GE.0)DD(KC2 )=0. 

IF (NO Cl SI (MCI ) .NE .0)DD( MCI ) =0 
IF< N00IS1 (MC2 > .NE.O )DD( MC2) =0 
WRITE (6, 168)1, 0D( MCI) ,DD(MC2) 

28 CONTINUE 

CALL POESS (OD ,1°,AL,NP,RAD,KC T ) 
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IF(KCT.FC.0)PEAD(5,100)ISLP 
IF(ISLP.EQ.0)G0T043 
K CK -G 

DG44 1=1, NP 
« C 1 = 2 «■ I P ( T ) - 1 

44 IF(DD( M C1 ) .LT.Ci )KCK = 1 

IF ( KCK .EQ .0 )G0T047 
WPI TE ( 6 , 1 75 ) 

IF(KCT,EG« n )REA0(5,104)C0F 

G0T043 

47 CALL FLIP (ISTR,DG,N0DIS1,KCT,C0F,ICY,IP,NP) 

IFIKCT.EQ.l ) G CTO 34 
IF ( IC V .E(J . 1 JG0T043 
I »R = 1 
ICY = 0 
G 0 T 0 4 6 

43 IF ( IC 1 .EQ .NOC )G0T01 5 

IC1 = IC1 ♦ 1 
9 E A D < 5 , 1 U 4 ) M 
WRITE (6,174 )W 
W =W*0 . 104 71976 
WRITE (6,152) 

GCT036 

16 IF ( IC .LT. TCASE ) GOTO 1 
1 f C FORK A T ( 12 ) 

1C1 P09KAT(4CH 
1 

1 C 2 r OR M A T ( 6 ; 4 ) 

1C2 FCRMAT(5tl2.« ) 

1 P 4 r ORKAT ( F 1 P • 0 ) 

105 FORMAT (6F10.0 ) 

IP 6 r 0RMAT(l4,4E12.4) 

107 FORMAT (F10»Q,E12.4) 

1C £ FORMA T ( 314 ,2F 10 .4 ) 

IF L FORMAT ( lh 1 ) 

151 FORMAT! 1H ) 

1 c - 2 FORMAT < 1 H l? ) 

1 F 3 FO p MAT (1H + ) 

154 FC d MAT(1X,’M=’,I4,5X,'N = ’,Ih,5X,'Ni«i=’,I < *,5X,'NE-',I4,5X, 

l'NC0N=',I4,5X,'NF0R = ’,I'») 

155 FCRMAT(1X,’PR= , ,L12.4,EX, , YM=’,E12.4,5X, , C0FX= , ,E12.4,5X, 
1'PH0= , ,E12.f,5X, , TEMP = , ,E12.4,5X, 'SPEED<RPM)=*,F10.2 ’ 

156 F0RMAT(1X, 'ELEMENT', 4X»’TYPE',7X, 'I', 7X,’J’,7X,*K’»7X,’L’ , 

110X, 'PI ' , 10X, 'RJ' ,7X, 'RKOBJ ' ,7X,'RL0 K ’ , 8X , ' BOBL ’ , 8 X , ' TE mP ' ) 

157 F0RMAT(6IS,6E 12.4) 

156 FORMAT ( 1 X , 'SPECIAL ELEMENTS */lH ,’ ELEMENT ', 10X , *PR ', 10X , ’Y h ' , 

1 8 X , 'COEX' , 9 X , 'RHO' ) 

155 F ORM A T ( 1 8 ,4E 1 2 • 4 ) 

16C FORM A T ( 1 X , ' 6 A SE MATERIAL PROPERTY VARIATIONS WITH TE MP ERA T LR E ’ ) 

161 FORMAT (6X, 'TEMP' » IPX » 'PR' ) 

162 FORMAT <6X ,' TEMP ', 10X, *Y* ' ) 

1 t 2 FORMAT <6X ,' TEMP ' ,8X ,' COEX ’ ) 

164 FORMAT (Fx TEMP ' ,9X ,* RHO ' ) 

165 FORM A T ( 1 X , ’ NODE ' , 6 X , * I II * , 6 X , ’ I V ' , 8 X , ’ III) ' , 8 X , ' V V ' ) 

166 F0RMAT(I5,2I8,2F10.4) 

167 FORMAT (IX, 'NODE', 10X,’rR’, 10X, ’FZ’) 

166 FORMAT* 15 ,2E12.4 ) 

166 FORMA T ( 1 H , ' N OD E ' , 1 0 X , ' UR ' , 1 C X , ' U 2 ’ ) 

17C FORMAT ( IX , 'ELMT ' ,4X , 'R STRESS’ ,4X,’T STRESS ', 3X ,' AX C TRESS',3X, 



1'SH STRESS*, 5X, *R , ,9X t f Z t ) 

171 FORMATUH , 'AVERAGE TANGENTIAL STRESS=' ,E12.4> 

172 r ORH A T { 1H 1 , 'NEW SPEED AND/OR FORCE S .W= * ,F1G *2 > 

1 ^2 FOPMATdhC, 'INTERFERENCE FORC^ VECTORS'/IH , • N ODE • , 6 X , • RAD IA L * , 
1 7 X , 'AXIAL * ) 

174 FORMAT ( lhl f 'SAME G E OM E T R Y . W = ' , F 10 . 2 ) 

175 FORMAT ( 1 HQ , 'LOSS OF RADIAL INTERFERENCE') 

176 FORMAT (IX , 'NODAL RESTRAINTS') 

1 7 7 roRMATdX, 'APPLIED NODAL FORCES') 

17e FORMATUX , 'NODAL DEFORMATIONS') 

179 FORMATdX , 'ELEMENT STRESSES') 

STOP 

END 


SUE ROUTINE AINVfRI,RJ,RK,RL,b,AI,AJ,BJ,BK,BL,DET,L,TRT,A) 

C INVERTS THE ARRAY A 

DIMENSION A ( 8 , 8 ) , TR I (200,10) 

CALL SFT21 ( A , 8 , 8 ) 

IF(TRI(L,2) .GE*O.C)GOTO 1 
ff I=TRI <L , 1 ) 

PJ= - T RI(L,2) 

RK=TRI (L , 7 ) 

R L = TR T ( L , 4 ) 

B =T R I < L , 5 ) 

B J = 0 * 

A I- 0.5*1 R I ♦ R K ) 

AJ = 0 • 5* ( R J + R L > 

A ( 2 , 1 ) - 1 / ( R I -R J ) 

A (6 ,5 ) “ M 2 , 1 ) 

A (2,2)- -M 2 , 1 ) 

A (6 ,6 ) - A (2,2 ) 

A <1,1 ) = Rj*A(2,?> 

A ( 5 , 5 ) = A ( 1 , 1 ) 

A < 1 ,2 ) = R I * A (2,1) 

4 (5,6)- A(l, 2) 

A (4,1)“ A(2,2)/8 
A (8,5)- A ( 4 , 1 ) 

A ( 4 , 2 ) = -A ( 4 , 1 ) 

A (8,6)” A ( 4 , 2 ) 

A (3,1 ) = R J* A ( 4 , ? ) 

A (7,5)- A ( 3 , 1 ) 

A (3,2)- R I* A (4,1) 

4(7,6)= A ( 3 , 2 ) 

A (4 ,3 )= 1/ (B* (RK-RL ) ) 

A ( S , 7 ) = A ( 4 , 3 ) 

A ( 4 , 4 ) = - A ( 4 , 3 ) 

A (8,8)= A ( 4 , 4 ) 

A (3,3 )= RL*A (4,4 } 

A (7,7)= A ( 3 , 3 ) 

A ( 3 ,4 ) = RK*A (4,3) 

A (7,8)= A ( 3 , 4 ) 

D ET = AI-AJ 
RETURN 

1 AI=TRI(L,1) 

AJ=TRI(L,?> 

DET=AT“AJ 
8 J= TR I (L , 7) 

BK= Tk I ( L , 4 ) 
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R L = TR I ( L , c ) 

9 - (6K+°L-BJ)/2 

D 1 - AJ-AI 

D2 = DIMtJ-BL) 

A ( 1 ,1 ) = A J/D 1 
A (5*5)- A ( 1 , 1 > 

A ( 1 ,2 ) = AI^BL/D' 5 
A (5 ,c ) = A ( 1 , 2 ) 

A < 1 ,4 ) = -AI*B J/D2 
A { 5 t 8 ) - A (1 ,4 ) 

A ( 2 ,1 ) = - 1/01 

A (6*5)- A(2,l ) 

A (2 , 2 ) = -FL/D2 
A (6,6)= A ( 2 , 2 ) 

A <2 ,4 ) = B J/D2 
A (6,6)= A (2 ,4 ) 

A (4,1)= 1/<6“*D1) 

A (8,5)= A (4 , 1 ) 

A (4,3)= -A ( h , 1 ) 

A (8 ,7 )= -A ( 4 , 1 ) 

A (3,1 )= -AJ*A (4,1) 
A ( 7 , E> ) = A ( 7 , 1 ) 

A ( 3 , 3 ) = - A ( 3 , 1 ) 

A (7,7)= - A ( 3 , 1 ) 

A (4 ,2 )= 1 /D2 
A ( 8 , 6 ) = A ( 4 , 2 ) 

A (4 ,4 ) = -A ( 4 , 2 ) 

A (8,6)= - A ( 4 , 2 ) 

A (3 ,4 ) = A I / D 2 
A (7,8)= A (7, 4) 

A (3,2)= - A ( 3 , 4 ) 

A ( 7 , 6 ) = -A (3, 4) 

RETURN 

END 


SUBROUTINE BCROUT(D,X,N,REPT,NW,NODlS) 

C ARROCEDURE TO APPLY CROUT REDUCTION TO A BANDED MATRIX 
C D WHICH HOLDS THE UPPER TRIANGLE OF THE SYMMETRICAL 
C BANDED MATRIX A SO THAT D(I , J - I ♦ I ) = A ( I , J ) . T h E LAST COLUMN 
C ( D ( I , J -I + 1 ) HOLDS the VECTOR B FOR THE SET OF EOLATIONS '.X=B* 

C ON EXIT THE VECTOR X HOLDS THE SOLUTION (X) TO THESE F DU « TION 5 . 
C TFE REDUCTION IS ARRANGED SO THAT IF CERTAIN SETS OF EQUATION' 

C ^ 1ST NOT BE SOLVED. (I. E. N0DIS(I) = 1) THEN THE APPROPRIATE RO*S 
C AND COLUMNS ARE NEGLECTED 

DIMENSION D (40C ,41) ,X (4 00 ) ,NODIS( 4G0 ) 

INTEGER REPT,CT,C1 
NW2=2*NW+ 1 
CT = C 

91 C T = C T + 1 

IF ( (NODIS(CT)-D.l) • G T . 0 ) GO TO 91 

if (rept.lt. i ) goto 25 

DO 92 1 = 1 , N W 2 

IF ( (NODIS (CT + I ) -0.1 ) . GT .0 ) G0 T 0 92 
D ( C T , T ♦ 1 ) = D(CT,I + 1)/D(CT,1) 

92 CONTINUE 
C i=CT+l 

no S 3 j r c 1 , h 
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IF { ( NODI S ( J )“ r, *l)*GT*0) GC TO 93 
SUM=Q. 

LI =J-1 
K 1 =J-NW2 

IF (Kl.LT.CT) Kl-CT 
DC94 K=K1 ? L1 

IF UNODISCK )-a.l).GT.O) GO TO 94 
SUM =SUM+ (D (K f J-K+l )**2)*D(K, J ) 

94 CONTINUE 

D ( J , 1 ) = D ( J f 1 ) “SUM 
IF (J.EO.N) GOTO 25 
L 2= J + 1 
N2 - J +N W 2 
IF ( N 2 * GE • N ) N2 = N 
DO 95 I “L 7 t N2 

IF ( < NODI S C I )-n»n,GT.O> GOTO 95 
SUM = 0. 

K 1 =I-NW2-1 
IF (Kl.LT.CT) K 1 = CT 
00 96 K = K 1 ,L1 

IF ( l I-K--.1 ) .GT.NW2) GOTO 96 

IF ( ( NODIS ( K ) -0 .1 ) .GT.O ) GOTO 96 

SUM =SUMMD(K v I-K*1)*D<K 9 J-K*1)*D(K 9 I) ) 

96 CONTINUE 

D ( J v I - J+ 1 ) = (D(JfI-J + l )-SUM)/0(J 9 l) 

95 CONTINUE 

93 CONTINUE 

25 NW3 - 2 *N W + 3 

D (CT,NW3)=D(CT f Nw3)/D(CT f l) 

C1=CT+1 
DO 97 I “ C 1 ,N 

IF ((NOriS(I) -0.1). GT.O) GOTO 97 
SUM =0. 

L1=I-1 
K 1 - I -N W2 

IF (Kl.LT.CT) K1=CT 
DO 98 K=K 1 , LI 

IF ( ( NOD I S ( K )-□.!) .GT.O) GOTO 98 
SUM -SUM+ (D(K 9 I-K + i ) *D ( K » 1 )*D(K*NW3) ) 

98 CONTINUE 

D ( I v N W 3 ) =(D(I f NU3)-SUM)/0(I t l) 

97 CONTINUE 

IF< (NODIS (N)-D* 1 ) .GT.O) GO TO 2 6 
X <N)=D(N ,NW3> 

26 L 1 - N - 1 

DO 99 1 = 1 t L 1 
L2=(N-I)+1 

sumzo.o 

L 3 = N -I 

IF ( (N0DIS(L3 ) -0 . 1 ) .GT.O) GOTO 99 
N 2 = L 3 +N W 2 

IF (N2.GT.N) N 2 =N 
DO 90 K=L2,N2 

IF ((NODIS (K >-0.1 ) .GT.O) GOTO 90 
SUM =SUM+(D(L3 f K-L3+l ) * X ( K ) ) 

90 CONTINUE 

X (L3)=0(L3>NW3)-$UM 

99 CONTINUE 
RETURN 
END 
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SUBROUTINE DIFDIS 

MKT, D,NODIS, MOVE, NOD, DD, DEL, ISTR, STR, NOOD»KCT> 

DIMENSION D(4GO,4l), N0DIS(400), MCVE(400>, 
100(400) , DEL ( 400 ) 

2 , IS TR (50,4) , S TR (50,2) 
COMMON/PLOCK2/MORO,NW221,NW22,NW 
CALL SETZ2 ( DE L , 400 ) 

CALL SETZ?(DD ,400) 

CALL SETZ4(NDD,400) 

CALL SE TZ 4 ( MO VE , 400 ) 

NW2 - MW22-1 
D 02 0 I = 1 » K T 

20 DD ( I ) = 0 ( I , NW?2+1 ) 

IF ( KCT.NE .0 ) G 0 T 090 
READ ( 5 , 1 ) NO CD 
WRITE (6,2 ) 

WRITE (6,251) 

0021 I = 1 ,NOCD 

READ(5, 252)0, K,IR, DR, IZ,DZ 

WRITE(6,252)J,K,IR,DR,IZ,DZ 

ISTR(I ,1 )=J 

I STR ( 1,2) =K 

I STR ( I * 3 ) - 1 R 

I STR ( 1,4 ) =IZ 

STR (1,1)- DR 

STR (I ,2)=DZ 

NDD(2*K-1) Z ( -2*0+ 1 ) *1 R 
NDD(2*K) = ( -2* J ) *1 Z 
DEL(2*K-1) = OR 

21 DEL ( 2 *K ) = DZ 

GOT 091 

90 D092 I- 1 , NOOD 

J -I STR ( I , 1 ) 

K =ISTR ( I ,2) 

DR=STR(I, 1) 

DZ-STP (1,2) 

IR=ISTR (1,3) 

IZ-ISTR ( I ,4 ) 

N DD ( 2 *K - 1 )=(-2*0+l)*IR 
NDD (2*K ) = (-2* J) * IZ 
DEL ( 2 *K -1 )-DR 
92 DEL ( 2 *K ) =QZ 

91 11=0 

D022 I = 1,KT 

MOVE(I) = I 

TF (NDDII ).LT.U) GOTO 70 
II = II + l 
MOVE ( I ) = II 
70 MCI = I-NW2 

MC2=I + NW2 

IF (MC1.6T.1) 01 = MCI 
IF (MC1.LE.1) 01 = 1 
IF (MC2.LT.KT) 02 = MC2 
IF (MC2.GE.KT) 02 = KT 
D 02 3 0 = 01,02 
I F ( NDD (O).GE.O) G0T023 

IF (O.GE.I) DD ( I ) = DD ( I ) -DEL (0 >*D(I ,0-1+1 > 

IF (O.LT.I) DO(I) = DD ( I ) -DEL (0 )*D ( 0 , I-J+l > 

23 CONTINUE 


N C D ( 4 0 0 > , 
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22 CONTINUE 

0024 I = 1,KT 
IF(NCDCI) .LT.C) GO TO 7 1 
IK = MOVE C I ) 

NFI = 0 
GCT07? 

71 IK = -NDD(I) 

NFI = 1 

IK = MOVE (IK) 

72 J1 = I 

MC2 = I+NW2 

IF (MC2.LT. KT) J2 = MC? 

IF (MC2.GE.KT ) JZ - K T 
DC 25 J : J1,J2 
IF(NDD( J) .LT.C) G0T07 4 
JK=MOVE ( J ) 

JJ = JK 

nfij = n 

NFJ = 0 
GCT073 

74 JK = -NDD(J) 

JK = MOVE (JK) 

NFJ = 1 

NFIJ = NF I*NF J 

73 MCI = NFJ+NFI-NFIJ 

IF (IK.GT. JKIDC JK,IK-JK+1 )=D(T, J-I+l ) +MC1*D( JK,IK-JK+1 ) 

IF ( IK .LE.JK)D(IK,JK-IK+1)=D(I,J-I+1)+MC1*D(IK,JK-IK+1 ) 

25 CONTINUE 
Jl-JJ+1 
J2-IK+NW2 

IF( J1.GT.J2)G0T075 
D 02 6 J = J1,J2 

26 D(IK,J-IK+1) =0 

75 CONTINUE 
DD(IK)=DD(I)+NFT*DD(IK) 

NODIS ( IK ) -NODIS(I )+NFI*NODIS ( TK ) 

24 CONTINUE 

I=MOVE(KT ) + l 
D 02 7 J=I,KT 
DC( J) -0. 

NODIS ( J>=0 
D 02 8 JJ-1 ,Nw22 

28 D (J,JJ)=0. 

27 CONTINUE 

D 02 9 I - 1 , K T 

29 D (I , N V ? ? 1 )-DD (I ) 

KT-MOVE(KT) 

WRITE (6,2 ) 

251 FORMAT! IX,’ DIFFERENTIAL DISPLACEMENTS ’/l ■ , ’NODE A’,5X,’N0CE B’,9 
1, ' I R ' ,°X, 'DR' ,9V, *12’ ,9X, * OZ ’ ) 

251 FORMAT ( IX ,It,5X,I6,5X ,16, ‘'l 1 • 4 , 5 X ,16, Ell. 4) 

1 FORMAT (515 ) 

2 FORMAT ( 1HH) 

RETURN 

END 
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6 F0PMAKT5 , 14X ,E12.4 ) 

7 FORMATUHD, f NON CON VE RGENCE f ) 
1C CONTINUE 

RETURN 

PhD 
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SUBROUTINE RECTfL,FW,NTRI,ALLK,U,TRI) 

DIMENSION NTRIl' , u0,4),T(S),FW(2 00) f Ul4C0),ALL«(4G0,41),F(o t p), 
1A (3 ,8 ) , TRK20G, ID) 

C OMMON/^LOCKi/TEMP, PR, YM, ALPHA, RHO 
IF(NTPI(L,4).E0.0)G0T08 

CALL A IN V < R I ,Rj,RK,RL ,B,AI,AJ,BJ,3K,3L,DET,L f TRI,A) 

G 0 T 0 v 

CALL T A I N V ( R I ,RJ,AI,AJ,BJ f BK f RL t DET,L,T c, I,A ) 

CONTINUE 

CALL SETZ1(E,8,4) 

IF ITPKL ,2).LT.0.0) GOTO 1 
IF (AI.Lt.O.C) GOTO 2 
ELF 3 A L 0 6 ( A J ✓ A I ) 

GOTO 3 
ELF = D 

PI 3 <<Ph * A j ) - (BL * AI) + ( R J * AT)) 

P 2 = -L-bK-BJ 
P 3 3 (RL-PK )/ (AJ-AI ) 

P4 3 RL-bK 
PS 3 AJ-AI 

81 = (P4**2-Ed**2)/(P5**2 ) 

B 2 = <P4**3-E J**3)/ (P5**3) 

P 3 = (3K**^-3* b K**2*AI*P3 + 3*b*<*<AI*P3)**2-R2*Ai**3> 

B 4 - (3*bK**2*P3-6* r 'K*AI*P3**?+3*'I**2* r 2) 

B 5 = {3*BK*P7**?-3*AI*B2> 

66 3 <BK**?-2*AI*FK*P3) 

El 3 f? + P 1 *ELF/°S 
F 2 - Pl+P^MA J+AI)*C.5 

E3 - P • c - * ( (B6+B1*AI**2)*£LF+2*BK*P4-2*AI*B1*P5 
1+C.5* C 1*(AJ**2-*I**2)) 

Eh - 0 *6* (BK**2*P5-2*AI*fcK*P4+feK*P4* ( A J + A I ) 

(P4*#2-tJ**2 )^ps/3 ) 

E 5 - f D ?*(AJ**3-AI**3)/3+Pl*(AJ**2-AI**2>/2>/P5 
E 6 3 P.5* (F6* (AJ**2-A I**2)/2+?*BK*P4* ( A J** 3- A 1**3 >/(3*F5) 

1 + Bl* (0.25*P5**4 + aI* <P5**3 )/3) ) 

E 7 3 (°2*0.25*(AJ**4-AI**4)+Pl*(AJ**3-AI**3)/3)/P5 
FP = <R3*ELF+?4*P5+B5*(Aj**2-Al**2)/2 

1 ♦ L 2 * < AJ**3-AI**3)/3)/3 

E 9 3 <R3*P5 + P4*< AJ**2 -AT**2 )/2 + 85*( A J**^-AI**3 )/3 
1 + 62* ( AJ**4-A1**4 )/4 ) / 3 
F 10 3 (~3*(AJ**?-AJ**2) /2+B4*(AJ**3-AI**3) /3 
1 + B5* < A J**4-AI**4 ) /4 + B 2* ( A J**5-A 1**5) /5 )/3 
Ell 3 <P2*(AJ**5-AI**5) /5 + Pl*( A J* * 4- A I **4 ) /4 ) /P5 
E 1 2 3 (P6*(AJ**3-AI**3)/3+2*&K*P4*<AJ**4-AI**4)/(4*P5) 

l + El*in.2*P5**5 + ri .5*AI*P6**4+AI**2*P5**3/3)>/2 
E 1 3 3 <P6*(AJ**4-AI**4)/4+2*BK*P4*(Ad**5-AI**5)/( c *PF) 

1 + 61* { (AJ**4*P5**2 ) /6-AI*(0.2* (AJ**5-AI**5) 
2-0.25*AT*(AJ**4-Al**4))/3))/Z 
Aul=( AJ-AI)*BK 
IFUJl.GT.O.r )G0T04 
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1 


5 


e> 


El=- El 
E2=~E2 
E 3- -E 3 
E 4 - -E 4 
Ebr-EB 
B6=-E6 
E 7- -E 7 
E8=-E8 
E 9- -E 9 
F 10--E 10 
E 11=-EU 
E12=-F12 
E 1 3 - -E 1 3 
GOTO 4 

p JR I - Rj-RI 
PLRK = RL-RK 
RLLN = ALCG(RL) 
p JL N = ALOG(RJ) 

RKRI r pk-PI 
RLRJ = RL-RJ 
RKLN r ALOG (RK ) 

RILN = ALOG(RI) 

E 2 = 8*(R JRI+ (RLRK-RJRI )/2> 
f-"4 = E**2* (RLRK+0.5*R JRI ) /3 

E 5 = H .5*B* (R J**2-~I**2+f J*RL D J-RI*RKPI+ (RLR j**?-PKRI**I ) / 3 ) 

F 6 = B**2*( (RJ**2-RI**2)/2 + 2*(RJ*RLRJ-RI* : -KRI)/3 
1+ <RLRJ**2-PKRI**2 ) /4 ) /2 
E 7 = n *(RJ**3-RI**3*1.5*(RJ**' , *RLRu-' r T**7*RKPI) 
l+RJ*PLPJ**2-RI*9KRI**2+(PLRJ**3-R^RI**3)/4)/o 
P9 = B**3* (RJRI+3-RLRK ) /12 

E 12 = B**?* (( Rj**3-RI ** 3 ) /2+RJ* *2 *ALiJ’R I** 2 *RKR1 
1+G.75*(RJ*RLRJ**2-RI*RKRI**2)+0.2*(PL P J**3-PKPI** 7 ))/’ 

IF < A p S ( R L R J ) •LT«l*£-08) GOTO 5 

r 1 - ( -PLRJ+RL*RLLN-R J*RJLN > / P LRJ 

E3-RLL*‘-C • 5 + ( Rj/kLRj) *{ 1-RJ* [ D LLN-KJLN )/RLRj ) 

F8-PLLN-i./3.MRJ/ p lRJ)*(0.5-(RJ/RLRJ>*(l-PJ*(RLL^-RJLN)/RLRJ)) 

E10=PJ**2/6 + RJ^ p LRJ*0.25 + RLRJ s!r *2*0.1 

Ell - (PL**5-RJ*»5 )*0 .05/9LRJ 

E13=RJ**4*D.5+4*kJ**3*RLRJ/3+1.5*Rj**' 5 *PlRJ** 2+C.°*RJ*RLRJ**3 
l+FLRJ**4/f 
GOTO 6 
El - P JL N 
E3-RJLN 
Efc = RJL^ 

E 10-R J**2 /6 
Ell = 0 • 2 5* R J ** 4 
E 1 3 = 0.S*RJ**4 

IF (A3SIRKRI) ,lT.1.E-Q3) GOTO 7 

El r R*{El-(-RKPI+RK*RKLN“RI*PlLN)/RKRI) 

E3-B**' 1 *{E3-RKLN + 0.5-(RI/RKRI)*(l-RI*(RKLN-RILN)/P^RI))/2 
EE=B**3*(E8-RKLN+l./3.-(RI/RKRl)*{0.5-(Pl/RKRI) 

1* (l-RI*(fcKLN-RILM)/RKRI ) ) )/3 
E1C = B**3*(E10-RI^*2/6“PI*RKR1*0 .25 -R.K R 1**2* 0 . 1 ) 

FU = B* (El 1- (RK**5-R 1**5 )*G.H5/RKRI ) 

E 1 3 -B **2 * (E13-RI**4*0.5-4*RI**3*RKftI/3-1.5*Rl**2* p K'I**2 
1-0.8*PI*RKRI**3-RKRI**4 /6 >*0.^5 
GOTO 4 
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DB<3,8>= R*DB(3*7) 

0B(4,3)r C1-2*PP)*0.5 

0B(4,4)r R*DB(4*3) 

0 B ( 4 T 6 ) “ DB<4,3> 

DB(4, 8) = DB(4,3)*Z 
IF< NTRI (L *4 ) • NE .0)G0T04 
D 03 1-1*4 
DB(I*4)rO. 

3 DB(I,8)=0. 

4 CONTINUE 

NR1 = NTRIIL, 1 ) *2 -1 
NR2 = NTRI<L,2>*2 -1 
NR3 = NTRI(L*3)*2 -X 
T D < 1 ) = U ( NR 1 ) 

T D ( 2 ) = U(NR2) 

TO ( 3 ) r u ( NR 3 ) 

T D < 5 ) - U (NR1 +1 ) 

T 0 ( 6 > = U ( NR 2 ♦ 1 ) 

T C ( 7 ) = U ( N R 3 + 1 ) 

IF(NTRI(L*4) .EG.0JG0T05 
NR4 = NTRI(L*4)*2 -1 
TD ( 4 ) = U (NR4 ) 

TDm = U ( N R4 ♦ 1 ) 

5 CONTINUE 

CALL STEP2(8,ST*T0TA t A,TD*DB, STRESS, AREA) 
WRITE (6*2 *0 )L , <$TRESS(I)*I=i,4) t R,Z 
2 5 C FORMAT ( 1 5 * 4E 12 . 4 * 2F 1 0 . 5 ) 

RETURN 

END 


SUBROUTINE SETZ1 ( ADUM *N *M ) 
DIMENSION A DU M ( N » M } 

DO 62 1-1 *N 
DO 62 Jri r M 
62 ADUH(IfJ) =0. 

RETURN 

END 


SUBROUTINE SETZ?(ADUM,N) 
DIMENSION AOUM(N) 

DO 62 1-1 * N 
62 ADUM C I ) -0. 

RETURN 

END 


SUBROUTINE SE TZ 7 ( N D UM , N , M ) 
DIMENSION NDUMtN,M) 

DO 62 ir 1 ,N 
DO 62 J=1 * M 
62 NDUM ( I * J ) ro 
RETURN 
END 


28 



SUBROUTINE SE TZ4 ( NDUM ,N ) 
DIMENSION NDUM(N) 

DO 62 T - 1 , N 
62 NDUM l I ) = P 
RETURN 
END 


SUBROUTINE SLIP(ISTR f DD,NODISl,KCT,COF,ICY f IP,NP> 
DIMENSION ISTR(50,4) r DD(4DO),N0DISl(4r0)fIP(50) 
IF(KCT.EQ.l ) G 0 1 0 1 1 
PE A D ( 5 , 1 ) COF 
11 ICY=ICY*1 

WRITE ( 6 , 5 ) I C Y 
WRITE (6,2 )C0F 
K CT = 0 

D 01 C 1 = 1, NP 
MC1=2*IP m-1 
MC2=MC1+1 

IF(N0DIS1(MC2).NE.0)G0T010 

RATI0=D01mC2>/DD(MC1) 

PATIO=ABS (RATIO) 

IF (RATIO.LT • C OF ) 60 T 0 1 0 
K CT - 1 

I STR ( I , 4 ) =0 
WRITE (6,i)IP(I) 

10 CONTINUE 

IF (KCT.EQ*Q)WRIT£(6 f 4 ) 

1 r OR M A T < F5 * 0 ) 

2 FORMAT ( 1H ^COEFFICIENT OF F kTC TI 0N= • , FS . 3 ) 

3 FORMATdh , f N0DE # ,I5 f * ALLOWED TO SLIP AXIALLY*) 

4 FORM A T ( 1 H f *N0 SLIPPAGE 1 ) 

5 FORMAT( 1N0, *SLI P CYCLE*,I3) 

RETURN 

END 


SUBROUTINE STEPl<MS,BTDB f A,ALLK t T t U,NTRI t L) 

C POSITIONS IN ALLK ARRAY 

DIMENSION BTDB(B f 8),A(8,8),ALLK(40Q,4l),T<P),U{4C,n),NTRI(2C0,4), 
1ELKA ( 8,8 ) 

C CM MON/PL OCK 1/TEMP, PR ,YM, AL PH? ,RHO 
M2 = MS/2 
M 3 = M2 

FACT = YM / ( ( 1 ♦ PR) * (1-2 * D P)) 

TEMP = YM*TEMP*A L PHA/ (1-2 *PR> 

D 01 I = 1,MS 
D 0 1 J = 1 t M S 
SUM = n.G 
D 02 K = 1 , M S 

2 SUM = SUM ♦ 6TDB <I,K)*A(K,J) 

ELKA(I,J) = SUM 
1 CONTINUE 

DG3 I = 1 * M S 
0 03 J = 1 ,MS 
SUM - n.G 
0C4 K = 1 ,MS 
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SUM = SUM + A (K ,1 )*ELKA (K , J ) 

BTDB(ItJ) = $UM*FACT 
CONTINUE 

IF (NTRI (L f 4) . EQ • G ) M 3 = 3 
DOS I = 1 , M 3 
1 1 = T ♦ u -M2 
M K - NT P I (L 1 1 I ) *2-1 
K K - MM ♦ 1 
II = I ♦ M2 

U(MM) = T (I ) * TE MP ♦ U ( M M ) 

U(KK) - T(II)*TEMP + U ( K h ) 

D 06 J = 1 , M 3 

JJ = J ♦ 4-M2 

NN = NTRI (L,JJ)*2-1 

LL = NN+1 

JJ z J + M2 

IF (MM.LE .NN) ALLK (MM,NN-KM+1) = A LLK < MM ,NN -MM+ 1 ) + BTDB<I f J> 

IF (MM .LE . LL ) ALLK (MM f LL-MM + 1 ) = AL LK ( MM , LL-MM + 1) ♦ BTDP(I,JJ) 

IF (KK ,LE .NN ) ALLK (KK,NN-KK+ 1 ) = A LLK ( K K , NN-KK + 1 ) ♦ BTDB ( 1 1 * J ) 

IF (HK.LE.LL) ALLK(KK,LL-KK + 1) = A LLK ( KK ,LL -KK + 1 ) ♦ BTDB(II,JJ) 

CONTINUE 

CONTINUE 

RETURN 

ENO 


SUBROUTINE STEP 2 { MS , ST ,TOTAfA,TD ,DB, STRESS, AREA ) 
DIMENSION A<8,8) f TD<8),DE(4,S>,STRESS(4> f AT0<8> 

C CMM ON /BLOCK 1 /TEMP, PR, YM, ALPHA, RHO 
FACT = YM/< <1+PP)*C1-2*PR>) 

D 0 1 1 = 1, MS 
SUM=C .0 
DC2 J=1,MS 

SUM= SUM + A ( I , J ) *TD ( J) 

AID < I > = SUM 
D 03 1=1,4 
SUM= C.G 
004 J = 1 , M S 

S UM = SUM + DB ( I , J )*ATC ( J ) 

STRESS (I )- FACT*SUM-TEMP*ALPHA*YM/ ( 1 -2*PR > 

STRESS(4)= F A CT*SU M 

ST= ST + STRE$S<?)*AREA 

T OT A = TOT ^ + AREA 

RETURN 

ENO 
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APPENDIX B 


A BRIEF CONSIDERATION OF THE STRUCTURE OF THE PROGRAM 

A flow chart of the computer program for a single case run is shown in figure 3. 

The functions of the program elements listed in appendix A will now be briefly described: 

AXIFE the main program element, which acts as executor, reading in data, 

calling appropriate subroutines, and routing according to options 
and multiple case processing requirements 

AINV inverts the shape function matrix (A) for trapezoidal elements 

BCROUT applies matrix reduction solution scheme to solve for deformations 

DIFDIS reduces system of equations according to interference (differential 

displacement) conditions 

PRESS calculates pressure distribution at interference location 

RECT prepares data for calculating element stiffness characteristics and 

thermal and inertial force effects 

RECTST prepares data for calculating element stresses 

SETZ1,SETZ2,1 

V sets elements of arrays of various dimensions and types to zero 

SETZ3,SETZ4 J 

SLIP allows axial slip of interfering nodes if appropriate 

STEP1 determines element stiffness characteristics and organizes them 

into the main stiffness matrix 

STEP2 calculates element stresses 

TAINV inverts the shape function matrix (A) for triangular elements 
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APPENDIX C 


SOME PROGRAM DETAILS 

Array dimensioning . - The following arrays would need redimensioning if the num- 
ber of nodes required exceeds 200: ALLK, TRI, NTRI, NODIS, FW, U, TU, DEL, DD, 
NDD, MOVE, ALLK1, NODIS1, and D. 

The following arrays would need redimensioning if the maximum nodal difference in 
an element (NW) were to exceed 19: ALLK and ALLK1 (see section Maximum nodal 
number difference) . 

The following arrays would need redimensioning if the number of constrained nodes 
were to exceed 50: NUV, ANUV, IP, AL, P, Z, B, and A. 

The appropriate modifications, should redimensioning be required, are apparent 
without further detailing. Other dimensioning limitations are unlikely to be significant. 
The storage requirement of the program, as listed, on a Univac 1100 machine is 53 004 
words. Note also that if redimensioning of arrays whose elements are set to zero is 
undertaken (through SETZ1 etc. ), then the call of the subroutines SETZ should be 
amended accordingly. 

Maximum nodal number difference . - The matrix reduction solution scheme after 
Crout (ref. 3) has been amended to a banded form for the purposes of the present pro- 
gram. The bandwidth of the solution is (2NW + 2), where NW is the maximum nodal num- 
ber difference in any element. In fact, the width of the main stiffness matrix ALLK, and 
its copy ALLK1, has been dimensioned as (2NW + 3), the nodal force vector being added 
on in the last column. 

Whilst it is apparent what the maximum nodal number difference is in an element 
where none of the nodes is ’’involved” in an interference, it is not so clear if this is not 
the case. If an interference is effected between components, then for elements involved 
in the interference the value of NW should be assessed as the maximum nodal number 
difference between nodes in the element and associated nodes in interference with the 
element. (Examples are shown in fig. 4. ) 

Equivalent loading . - If a continuously distributed load is applied to a structure, the 
loading must be reduced to an equivalent nodal distribution. An example is shown in 
figure 5. 

Frictional condition at interference surfaces . - If two nodes have an imposed radial 
interference, they may also be constrained axially or, alternatively, allowed to move 
freely axially (to slip). If no slip is permitted between two radially constrained nodes, 
then equal and opposite forces will be set up on the nodes, there being an axial compo- 
nent. The reality of whether or not axial slip occurs between the nodes might be thought 
of as whether the ratio of axial to radial force exceeds some assumed coefficient of 
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friction. An arrangement for automatically allowing slip has been incorporated into the 
program (cards T and U). 

The option to allow slip to occur should only be chosen when all nodes in interference 
are located at the same radius and are axially adjacent. With the slip option effective, 
all the node pairs should initially be constrained axially. If slip of any node pair actually 
occurs, then a recalculation of ciistortioos, etc. , is carried out. Several cycles may be 
required until either no slip occurs in the remaining nodes constrained axially or all 
nodes have been freed to slip axially. 

The process of slip can be investigated by manual intervention and, although this is 
more time consuming, it jnay be mpre appropriate in some circumstances. 
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APPENDIX D 


AN EXAMPLE 


To demonstrate the use of the program, an example will now be presented. This is 
the inner ring of a roller bearing having an interference fit with its hollow shaft. An ex- 
tremely course mesh of elements has been chosen (see fig. 6) to avoid too much detail. 
Figure 6 shows node and element numbering chosen together with dimensions given in 
meters. A symmetrical situation is considered and therefore only half the geometry has 
been considered. The following data were chosen for the example: 


Rotational speed, rpm 

Ring and shaft material 

Poisson's ratio . 

2 

Young’s modulus, N/m 

Coefficient of expansion, °C~* 

Density, kg/m 

Ring temperature (relative to the shaft), °C 
Shaft temperature (no thermal expansion), °C 

Radial interference, m 

Coefficient of friction (for axial slip) .... 


. 5000 
. steel 
. . 0.3 
2X10 11 
. IQ" 5 
. 7800 
. . 20 
. . . 0 
4xl0" 5 
. 0.5 


Note that nodes 1 to 7 are restrained axially to fulfill symmetry requirements and that 
the nodes in radial interference 4/5, 11/12, 20/21, 29/30, and 38/39 are at the same 
radial location; that is, there are two nodes at the same spatial position. 

It should be noted that the technical aspects of this example are not the concern of 
this report. The example is solely to demonstrate a typical data input and the form of 
the output. The example used 13 seconds of central processing unit time on a Univac 
1100 computer. The data input and program output are now presented. 


35 



1 







CAPO 

A 

ROLLER BEARING INNER 

KING/SI UNITS 



CARO 

B 

1 







CARP 

C 

28 

47 11 

7 





C A 9Q 

D 

. 

300 OE GP 

• 2 00 0 1 

12 

• 1000E-04 • 78 0 0 E 

04 


r s n n 

E 


5000 . 






CARD 

F 

1 

1 «> 

8 2 

1 




CARD 

G 


.064 

.067 


.003 

.003 

2 0 . 

C A R 0 

H 

7 

1 1 1 

IP 4 

3 







.05 8 

.061 


.003 

.003 

20 . 



4 

1 1 3 

12 6 

5 







.053 

.058 


.003 

.003 




5 

1 14 

1 3 7 

6 







.048 

.053 


.003 

.003 




f 

1 IP 

1 7 9 

8 







.064 

.067 


.003 

.003 

20. 



p 

1 20 

19 11 

10 







.05 8 

.06 1 


.003 

.003 

?C . 



9 

1 22 

21 13 

12 







.053 

.058 


.003 

.00 3 




1 0 

1 23 

22 14 

13 







.048 

.053 


.003 

.003 




1 1 

1 25 

24 16 

15 







• 06R 

.071 


• 003 

.003 

2° . 



1 2 

1 26 

25 17 

16 







.067 

.069 


. 00 3 

.003 

20. 



1 7 

1 27 

26 18 

17 







.064 

.067 


.003 

.003 

20. 



1 5 

1 29 

28 20 

19 







.05 8 

.061 


.003 

.003 

20. 



1 6 

1 31 

30 22 

21 







.053 

.058 


.003 

.003 




I 7 

1 37 

31 23 

22 







. 04 e 

.05 3 


.003 

.003 




1 fi 

1 34 

33 25 

24 







.06 9 

.071 


.003 

.003 

20 . 



1 9 

1 35 

34 2b 

25 







.067 

.069 


.003 

.003 

20 . 



2 C 

1 36 

35 27 

26 







.064 

.067 


.003 

.003 

20 . 



22 

1 3 P 

37 29 

28 







.05 8 

.061 


.003 

.003 

2 0. 



23 

1 40 

39 3 1 

30 







. 053 

.058 


.003 

.003 




2 4 

1 41 

40 32 

31 







.048 

.053 


. DC 3 

.003 




25 

1 43 

4 2 40 

39 







.053 

.058 


.003 

.003 




26 

1 44 

43 4 1 

40 







. 048 

.053 


.003 

.003 




27 

1 46 

45 43 

42 







.053 

.058 


.00 3 

.003 




2 8 

1 47 

46 44 

43 







. C4e 

.05 3 


.00 3 

.003 
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CAPD L 


1 C 1 

? C 1 

3 0 1 

4 0 1 

5 0 1 


t> u 

1 











7 C 

c 

1 











4 



5 


1 

- . 4 COPE 

-0 4 


1 



11 



12 


1 

- • 4 COPE 

-0 4 


1 



2D 



21 


1 

- • 4 000 E 

-0 4 


X 



2® 



30 


1 

- • 4 CODE 

-04 


1 



3 fl 

l 



7 9 


1 

- • 4 GO 0 E 

-04 


1 



X 

5 

5 8C PE 

-0 1 










5 

300 OE 

-02 










1 2 

3C00E 

-02 










21 

3 00 PE 

-02 










3D 

3C00E 

-02 










39 

• 6 

3C00E 

-0 2 










ROLLER 

BEARING INNER RING/SI UNITS 






Mr 9 3 

N - 

47 

NW r 

1 1 

NE r C 

NC ON r 

7 

NFO R- 0 


PR r .3000 + 00 

YM: 

: .2000+12 

COEXr 

. 1000-04 

RHO= .7800+0M TEMP 

r .0000 SPtEO < R P M ) 

lLEMENT 

‘ TYPE 

I 

J 

K 

L 

R I 

RJ 

RKOFJ 


RLOBK 

BCBL 

1 

1 

9 

8 

2 

1 

.6400-01 

.6700-01 

.0000 


.3000-02 

.3000-02 

2 

1 

10 

9 

j 

2 

.6100-01 

.6400-01 

.0000 


.3000-02 

.3000-02 

3 

1 

11 

10 

4 

3 

.5800-01 

.6100-01 

.0000 


.3000-02 

.3000-02 

4 

1 

13 

12 

6 

5 

.5300-01 

.5800-01 

.0000 


.3000-02 

.3000-02 

S 

1 

14 

13 

7 

6 

.4800-01 

.5300-01 

.0000 


.3000-02 

.3000-02 

6 

1 

18 

17 

9 

3 

.6400-01 

.6700-01 

.0000 


.3000-02 

.3000-02 

7 

1 

19 

18 

10 

9 

.6100-01 

.6400-01 

.0000 


• 30 CQ - 02 

.3000-02 

9 

1 

20 

19 

11 

10 

• c 8 00 -0 1 

.6100-01 

.0000 


.3000-02 

.3000-02 

o 

1 

22 

21 

13 

12 

.5300-01 

.5800-01 

.0000 


.3000-02 

.3000-02 

10 

1 

23 

22 

14 

13 

.4800-01 

.5300-01 

.0000 


.3000-02 

.3000-02 

11 

1 

25 

24 

16 

15 

.6900-01 

.7100-01 

.0000 


.3000-02 

.3000-02 

12 

1 

26 

25 

17 

16 

.6700-01 

.6900-01 

.0000 


.3000-02 

.3000-02 

1 3 

1 

27 

26 

18 

1 7 

.6400-01 

.6700-01 

.0000 


.3000-02 

.3000-02 

14 

1 

28 

27 

19 

18 

.6100-01 

.6400-01 

.^000 


.3000-02 

.3000-02 

15 

1 

29 

28 

20 

19 

.5800-01 

.6100-01 

.0000 


.3000-02 

.3000-02 

16 

1 

31 

30 

22 

21 

.5300-01 

.5800-01 

. 0000 


.3000-02 

.3000-02 

17 

1 

32 

31 

23 

22 

.4600-01 

.5300-01 

."000 


.3000-02 

.3000-02 

18 

1 

34 

33 

25 

24 

.6900-01 

.7100-01 

• 0000 


.3000-02 

.3000-02 

19 

1 

35 

3M 

26 

25 

.6700-01 

.6900-01 

.0000 


.3000-02 

.3000 -02 

20 

1 

36 

35 

27 

26 

.6400 -01 

.6700-01 

.0000 


.3000-02 

.3000-02 

21 

1 

37 

36 

28 

27 

.6100-01 

.6400-01 

.0000 


.3000-02 

.300Q-U2 

22 

1 

38 

37 

29 

28 

.5800-01 

.6100-01 

.0000 


.3000-02 

.3000-02 

23 

1 

MO 

39 

31 

30 

.*300-01 

.5800-01 

.0000 


.3000-02 

.3000-02 

24 

1 

Ml 

40 

32 

31 

.4800-01 

.5300-01 

.CCOO 


.3000-02 

.3000-02 

25 

1 

M3 

42 

40 

39 

• c 3 00 -01 

.5800-01 

.0000 


.3000-02 

.3000-02 

26 

1 

MM 

43 

41 

40 

.4800-01 

.5300-01 

.0000 


.3000-02 

.3000-02 

27 

1 

M6 

45 

43 

4? 

.5300-01 

.5800-01 

.0000 


.3000-02 

.3000-02 

28 

1 

M7 

46 

44 

43 

.*800-01 

.5300-01 

.0000 


.3000-02 

.3000-02 


CASH N 
CAPO 0 
CAPD P 


CAPD 0 
CAPD R 
CA^D S 


CARD T 
CARD U 


5COC.CC 


TEMF 


• 2000 + D 2 

• 2 C 0 0 +0 2 
. 2GG0 +02 
. 0000 
.GCOO 
. 2000+02 
. 2 C DC + 02 

• 20 CO +02 

• COCO 

• CCCO 

• 2CG0 +02 
. 2000+02 
. 2000+02 
• 2000+02 
• 2000+02 

• CCOO 

• COOO 
.2000 +02 
. 2000+02 
. 2000+02 
. 2000+02 
. 2000+02 
.CCOO 

• CCOO 

• C Q 0 0 

.COOO 

• COOO 

.r ooo 
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NODAL 

RESTRA IMS 




NODL 

IU 

IV 

uu 

vv 

1 

C 

1 

• CD00 

• GOOD 

2 

0 

1 

.0000 

. 0000 

3 

0 

1 

.0000 

• GO 00 

4 

3 

1 

.0000 

.0000 

5 

0 

1 

.0000 

.0000 

o 

0 

1 

.0000 

. 0000 

7 

0 

1 

.0000 

.0000 


applied nodal forces 

NODE FR FZ 

DIFFERENTIAL DISPLACEMENTS 


NOD E A 

NODE B 

IR 

D ° 

IZ 


4 

5 

1 

-.4000-04 

1 

.0000 

11 

12 

1 

-.4000-04 

1 

.0000 

20 

2 1 

1 

-.4000-04 

1 

.0000 

29 

30 

1 

-.4000-04 

1 

.0000 

36 

39 

1 

-.4000-04 

1 

.0000 


NCDAL 

DEFORMATIONS 


NODL 

UR 

UZ 

1 

.3274-04 

. COOO 


.3241-04 

.0000 

3 

.321 1 -04 

.CCOC 

4 

.3185-04 

.COOO 

5 

-.8153-05 

.COOO 

o 

-.8221-05 

. COOC 

7 

-.6386-05 

. COOO 

8 

.3273-04 

- . 3365-06 

9 

.3240-04 

- . 298C-06 

ID 

.3210-04 

- .2587-06 

11 

.3184-04 

- . 226 9-06 

12 

-.8156-05 

- .2269-06 

13 

-.8221-05 

- • 208 9-06 

14 

-.6387-05 

-.2185-06 

15 

.3311-04 

- . 0061-06 

16 

.3288-04 

-.7463-06 

1 7 

.3266-04 

-.6843-06 

1 

.3236-04 

— . 605 8 — 06 

19 

.3209-04 

- . 5238-06 

20 

.3184-04 

- . 4509-06 

21 

-.8160-05 

-.4509-06 

22 

-.8219-05 

- .4155-06 

23 

-.8379-05 

- . 4412-Qo 

24 

.3302-04 

- . 1169-05 

28 

.3278-04 

1106-05 

2 6 

.3257-D4 

- . 1041-05 

27 

.3729-04 

- . r 42 8-C 6 

28 

.3207-04 

- . 308C-0fc 

29 

.3184-04 

-.6639-C6 

30 

-.8163-05 

-.6639-06 
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31 

-.*196-05 

-.6215-06 

32 

-.8350-05 

- .6669-06 

33 

• 3 2°2 -04 

-.1533-05 

34 

.3269-04 

-.1466-05 

35 

.3247-04 

-.1402-05 

36 

.3219-04 

-.1295-05 

37 

.3197-04 

- . 1166-05 

3b 

.3189-04 

-.6572-06 

39 

-.8113-05 

-.6572-06 

40 

-.8134-05 

- .5265-06 

41 

-.8291-05 

- . 683 3-06 

42 

-.7822-05 

-.9279-06 

43 

-.8029-05 

1019-05 

44 

-.8220-05 

-.1073-05 

4 5 

-.7714-05 

-.1036-05 

46 

-.7926-05 

1159-05 

47 

-.8161-06 

-.1246-05 


tLEMENT STRESSES 


ELM T 

R STRESS 

T SThESS 

AX STRESS 

SH STRESS 

R 

Z 

1 

-.9436+06 

.S87C+C8 

-.1521*07 

-.2004*06 

.06550 

.00150 

2 

-.3029+07 

.6110*08 

-.4024+07 

-.2879*06 

.06250 

.00150 

3 

-.5279+07 

.6404*08 

-.6183*07 

-.3000+06 

.05950 

.00150 

4 

-.5069+07 

-.2981*08 

.4065*07 

-.1027*06 

.05550 

.00150 

5 

-.1803+07 

-.3222*08 

.4039*07 

.6720*05 

.05050 

.00150 

6 

-.1724+07 

.5850*08 

-.1113+07 

-.2380*06 

.06550 

.00150 

7 

-.3869+07 

.6086*08 

-.3805+07 

-.8573*b6 

.06250 

.00150 

6 

-.5967+07 

.6376*08 

-.6361*07 

-.1067*07 

.05950 

.00150 

9 

-.5409+07 

-.3001*06 

.3726*07 

-.3891+06 

.05550 

.00150 

1C 

- • 1 94 1 +D7 

-.3224*08 

.4058*07 

.1360+06 

.05050 

.00150 

11 

-.6136+06 

.5398*08 

.9208+05 

.1390*05 

.07000 

.00150 

12 

-.1971+07 

.5566*08 

-.3603+04 

-.2875+05 

.06800 

.00150 

13 

-.3296+07 

. 5803*08 

-.4447+06 

-.1742+06 

.06550 

.00150 

14 

-.6111+07 

. 6030*08 

-.3037*07 

-.1631+07 

.062S0 

.00150 

15 

-.7209+07 

.6327*06 

- .6606 + 07 

-.2492*07 

.05950 

.00150 

16 

-.6414+07 

- . 305 6 + 08 

.2876*07 

-.8247+06 

. G5 5 6 0 

.00150 

17 

-.2253+07 

-.3226*08 

.4029+07 

-.8681+05 

.05050 

.00150 

Id 

-.6178+06 

.5370*08 

.7989+05 

-.1152+05 

.07000 

.00150 

19 

-.1928+07 

.5542*08 

.9309+05 

.3017+05 

.06800 

.00150 

20 

-.3920+07 

.5766*08 

-.1192*06' 

-.4837+05 

.06550 

.00150 

21 

-.7082+07 

• 606 1 + 08 

-.2661*06 

-.7896+06 

.06260 

.00150 

22 

- .1355+08 

. 6109+08 

-.7345+07 

-.5204+07 

.05950 

.00150 

23 

-.7756+07 

- .3123*0e 

.1577+07 

-.2016+07 

.05560 

.00150 

24 

-.2440+07 

-.3229*08 

.3626+07 

-.7835+06 

.05050 

.00150 

25 

-.5458+07 

-.3122+06 

-.2213+07 

-.4622+0? 

.05550 

.00150 

26 

-.1945+07 

-.3218+08 

(.2520+07 

-.1416+07 

.06050 

.00150 

27 

-.3067+06 

- .2858 + 08 

-.3939+06 

-.1048+07 

.05550 

.00150 

26 

-.9778+06 

- . 321 6*08 

.4808+06 

-.9920*06 

.05050 

.00150 

AVERAGE 

TANGENTIAL 

STRESS: 

.7186*07 
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INTERFERENCE FORCE VECTORS 


NODE 

RADIAL 

A XI AL 

5 

.5674+03 

. COOO 

12 

.1175+04 

.9051+02 

21 

.1296+04 

.2456+03 

30 

.1639+04 

.4642+03 

39 

.2036+04 

. 1349 + 04 

NODE 

12 

21 

30 

39 

INTERFERENCE PRESSURE (RADI ALLY .R^CONST) 
.6484+07 
.6634+07 
.7755+07 
.6419+07 
.2907+08 


SLIP CYCLE 1 

COEFFICIENT OF FRICTICN= .500 
NODE 39 ALLOWED TO SLIP AXIALLY 


INTERFERENCE FORCE VECTORS 
NODE "ADIAL AXIAL 

5 .5590+03 .COOO 

12 .1160+04 .1254+03 

?1 .1304+04 .3115+03 

30 .1794+04 .1671+04 

39 .1910+04 .COOO 

NODE INTERFERENCE 

t 

1 2 
21 
30 
39 

SLIP CYCLE 2 
COEFFICIENT OF FRICTICN= .500 


NODE 

30 ALLOWED 

TO SLIP AXIALLY 

INTERFERENCE FORCE 

VECTORS 

NODE 

RADIAL 

AXIAL 

5 

.5514+03 

. COOO 

12 

.1146+04 

. 2261+03 

21 

.1 36 7+04 

. 1740+04 

30 

.1859+04 

.COOO 

39 

.1820+04 

.ccoc 


40 


PRESSURE (RADIALLY. RrCO NS T) 

.6379+07 

.6564+07 

•756C+07 

.8053+07 

.2659+08 



NODE INTERFERENCE PRESSURE (RADIALLY .RrCONST ) 


5 


•6312+07 

12 


.6418+07 

21 


.7943+07 

30 


.8756+07 

39 


.2497+08 

SLIP 

CYCLE 

3 

COEFFICIENT 

OF FRICTION- .500 

NODE 

21 ALLOWED TO SLIP AXIALLY 


INTERFERENCE FORCE 

VECTORS 

NODE 

RADIAL 

AXIAL 

5 

.5436+03 

.COOO 

12 

.1196+04 

. 1624 + 04 

21 

.1411+04 

.COOO 

30 

.1831+04 

.COOO 

39 

.1776+04 

.COOO 


NODE 

INTERFERENCE PR ES SURE ( R ADI ALL Y . R-CONS T ) 

b 

.6066+07 

12 

.6793+07 

21 

.8249+07 

30 

.8593+07 

39 

.2436+08 

SLIP 

CYCLE 4 

COEFFICIENT OF FRICTION- .500 

NODE 

12 ALLOWED TO SLIP AXIALLY 


INTERFERENCE FORCE 

VECTORS 

NODE 

RADIAL 

AXIAL 

5 

.5921+03 

. COOO 

12 

.1229+04 

.0000 

21 

.1385+04 

.COOO 

30 

.1811 +04 

.COOO 

39 

.1754+04 

.COOO 

NODE 

INTERFERENCE PR ES SURE ( R «D I A LL Y . R = CONS T ) 

b 


.6756+07 

1 2 


.6955+07 

21 


.8030+07 

?U 


.8536+07 

39 


.2403+08 


SLIP CYCLE 5 

COEFFICIENT OF FRICT I CN- .500 
NO SLIPPAGE 
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NODAL 

DEFORMATIONS 


NODE 

UR 

UZ 

1 

.3236-04 

.cooc 

2 

.3201-04 

.cooo 

3 

.3173-04 

• cooo 

4 

.3152-04 

.cooo 

5 

-.8479-05 

.cooo 

6 

-.8640-05 

.cooo 

7 

-.8801-05 

.cooo 

6 

.3239-04 

-.2722-06 

9 

.3204-04 

-.3145-06 

10 

.3176-04 

-.3577-06 

11 

.3156-04 

-.4189-06 

12 

-.8442-05 

- .344 1-07 

13 

-.8596-05 

- . 1693-06 

14 

-.8754-05 

-.3182-06 

15 

.3290-04 

-.4061-06 

16 

.3267-04 

-.4749-06 

17 

.3245-04 

- . 5589-06 

18 

.3213-04 

- .6357-06 

19 

.3186-04 

-.7110-06 

20 

.3167-04 

- . 8327-06 

21 

-.8331-05 

-.7011-07 

22 

-.8464-05 

-.3411-06 

23 

-.8609-05 

- .6330-06 

24 

.3300-04 

-.7711-06 

25 

.3277-04 

- .8355-06 

26 

.3255-04 

-.8945-06 

27 

.3227-04 

-.9702-06 

28 

.3203-04 

-.1057-05 

29 

.3185-04 

-.1237-05 

30 

-.8146-05 

-.1114-06 

31 

-.8228-05 

-.5216-06 

32 

-.8361-05 

-.9313-06 

33 

.3309-04 

-.1136-05 

34 

.3286-04 

-.1193-05 

35 

.3265-04 

-.1250-05 

36 

.3239-04 

-.1327-05 

37 

.3222-04 

-.1426-05 

38 

.3219-04 

-.1637-05 

39 

-.7808-05 

-.1700-06 

40 

-.7878-05 

- . 7062-06 

41 

-.8031-05 

-.1188-05 

42 

-.7265-05 

-.2551-06 

4 3 

-.7466-05 

-.8584-06 

44 

-.7657-05 

1394-05 

45 

-.6871-05 

-.3558-06 

46 

-.7071-05 

-.9782-06 

47 

-.7292-05 

-.1554-05 



ELEMENT 
ELM T 

STRESSES 
R STRESS 

T STRESS 

AX STRESS 

SH STRESS 


Z 

1 

-•8323+06 

. 5700 + 08 

-.3595+07 

-.2358+06 

.06550 

.00150 

2 

-.2811+07 

.6117+08 

-.8954+05 

-.2859+06 

.06250 

.00150 

3 

-.5270+07 

. 660 1 + 08 

.4106+07 

-.1301+06 

.05950 

.00150 

4 

-.5352+07 

- .3387 + 08 

-.4976+07 

.1540+05 

.05550 

.00150 

5 

-.1917+07 

-.3332+08 

.5679+07 

-.6740+04 

.05050 

.00150 

6 

-.1602+07 

.5708+08 

-.3089+07 

- .4264 + 06 

.06550 

.00150 

7 

-.3335+07 

.6120+08 

-.1539+06 

-.9331+06 

.06250 

.00150 

& 

- .5822+07 

.6593+08 

.3608+07 

-.3761+06 

.05950 

.00150 

V 

-.5866+07 

-.3373+08 

-.4963+07 

.3477+04 

.05550 

.00150 

10 

-.2162+07 

-.3303+08 

.5662+07 

-.1655+06 

.05050 

.00150 

11 

-.5979+06 

. 536 7 + 08 

.1035+06 

-.1777+05 

.07000 

.00150 

12 

-.1992+07 

. 5506 + 08 

-.8765+06 

.1610+06 

.06800 

.00150 

13 

-.2902+07 

. 575 3 + 06 

-.1278+07 

-.1042+07 

.06550 

.00150 

14 

-.5091+07 

. 6094 + 08 

-•5596 + 0t> 

- .1778 + 07 

.06250 

.00150 

15 

-.7086+07 

.6573+08 

.2596+07 

-.6713+06 

.05950 

.00150 

16 

-.7185+07 

- . '*346 + 08 

-.4803+07 

-.1672+06 

.05550 

.00150 

1 7 

-.2525+07 

-. 3245+06 

.5469+07 

-.8044+06 

.05050 

.00150 

1 6 

-.6332+06 

. 5391 + 08 

.3459+05 

.1969+05 

.07000 

.00150 

1 9 

-.1903+07 

. 5559 + 08 

-.1542+06 

-.1607+06 

.06800 

.00150 

20 

-.4396+07 

.5773+08 

-.2527+06 

-.8113+06 

.06550 

.00150 

21 

-.8618+07 

.6046+08 

-.2579+06 

-.1600+07 

.06250 

.00150 

72 

-.1374+08 

6389+08 

.6537+06 

-.1744+07 

.05950 

.00150 

23 

-.791 3+07 

- . 3246 + 06 

-.4009+07 

-.1535+07 

.05550 

.00150 

24 

-.2350+07 

-.3152+08 

.4554+07 

-.1868+07 

.05050 

.00150 

25 

-.3973+07 

-.2921+08 

-.2042+07 

-.3472+07 

.05550 

.00150 

26 

-.1557+07 

- .3049 + 08 

.2334+07 

-.2233+07 

.05050 

.00150 

27 

.1353+06 

-.2590+08 

- .3783 + 06 

-.7005+06 

.05550 

.00150 

23 

-.3622+06 

-.2917+08 

.4595+06 

-.1197+07 

.05050 

.00150 

AVERAGE 

TANGE NT I AL 

STRESS- 

. 7186+07 
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TABLE I. - DATA CARD SEQUENCE FOR A SINGLE CASE RUN 


Card(s) 

Comment 

A 

= 1 (number of geometries) 

B 


C 

= 1 (number of cases for the geometry) 

D 


E 


F 


G,H 

repeat as appropriate 

I 

if there are "special” material elements, repeat as appropriate 

J,K 

if base material properties are temperature dependent, repeat as appropriate 

L 

if nodes are constrained, repeat as appropriate 

M 

if nodal forces are applied, repeat as appropriate 

N 

if = 0, terminate data; if = 1, continue for an interference case 

o 


P 

repeat as appropriate 

Q 

if = 0, go to card T 

R 


S 

repeat as appropriate 

T 

if = 0, terminate data; if = 1, continue for slip analysis 

u 





Type 1 Type 2 




Triangular axisymmetric element 

Figure 1. - The three types of element incorporated into the computer program. 
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Radial (r) 


I ► Axial (z) 

Figure 2. - The coordinate system. 



Figure 3. - Flow chart for a single case run. 
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Element not Involved in an Elements having interference at Cylinder with uniformly distributed compressive loading of Q/unit 

interference. MW ■ nodes 2-3 and 16-17. Element (T) length on |t s ou t er surface 

16-1-15 NW -17-1-16; element (T) 

NW - 18 - 2 - 16 Equivalent forces at - 

Nodes 1 and 4 : - 

Figure 4. - Determination of nodal number difference (NW) in element. /7 tKL0 L 1\ 3ttRL0 

Nodes 2 and (— + Q ^trR “ 2 * 2 /“ — 4 — 

Figure 5. - Determination of node applied forces due to distributed load. 



Figure 6. - The example: the inner ring and shaft of roller bearing. (Not to scale. ) 
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